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Abstract 

Techniques for using Gaussian QCD sum-rules to predict hadronic resonance properties are developed for 
single-resonance and two-resonance phenomenological models, and criteria are developed for determining which 
of these models is required for analyzing a particular hadronic channel. The vector current sum-rule coupled 
to the p meson is shown to be consistent with a single resonance model, and the Gaussian sum-rule analysis 
results in an accurate p mass prediction which exhibits excellent agreement between the theoretical prediction 
of the Gaussian sum-rule and the phenomenological model. A two-resonance model is shown to be necessary for 
the Gaussian sum-rule for the non-strange quark scalar (fin) currents. The two-resonance Gaussian sum-rule 
analysis of the isoscalar and isovector (/ = 0, 1) fin scalar mesons exhibits excellent agreement between the 
theoretical prediction and phenomenological model. The prediction of the resonance properties of the / = 0, 1 
fin scalar mesons in this two-resonance model provides valuable information for the interpretation of the scalar 
mesons, including the A(1775). 



1 Introduction 

QCD sum-rules are an established technique for relating hadronic properties to theoretical QCD predictions. The 
most frequently used sum-rules for this purpose are the Laplace |jj and finite-energy sum-rules j^j . The significance 
of using the finite-energy sum-rule (FESR) as a supplementary constraint on a Laplace sum-rule analysis was most 
clearly established through the Gaussian sum-rules || 

G(s,r) = iJ^=e X J-^L) P (t)dt , r>0 (1) 

to ^ ' 

where pit) is a hadronic spectral function with physical threshold to- The quantity G (s, t) on the left-hand side of 
is determined from a theoretical calculation of a correlation function II (Q 2 )- The composite operators used to 
construct the correlation function serve as interpolating fields for the hadronic channel represented by pit) on the 
right-hand side of (Q). 
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The Gaussian kernel in ([!]) implies that G (s, r) satisfies a diffusion equation 

d 2 G(s,r) _ 8G(s,t) 
ds 2 dr 



(2) 



The relation between Gaussian and finite-energy sum-rules was established m by showing that the resonance plus 
continuum model for p(t), when evolved through the diffusion equation, would only reproduce the QCD prediction 
at large energies (large r) if the resonance and continuum threshold sq were related by the n = member of the 
following FESR family. 

F n (so) = lJt n P (t)dt , n = 0,1,2,... (3) 

to 

As with the Gaussian sum- rule ([!]), the left-hand side of ^ is obtained from a QCD prediction. 

Sum-rules with Gaussian-like (Lorentzian) kernels have already proven to have strong predictive power in nuclear 
physics Q . The purpose of this paper is to demonstrate that Gaussian sum-rules can also be used predictively in 
hadronic physics, and to develop techniques for the analysis of Gaussian sum-rules in hadronic physics. Section |^ 
reviews the formulation and conceptual advantages of Gaussian sum-rules, and develops techniques for analyzing 
a single-resonance plus continuum phenomenological model. Section || applies these techniques to the p meson 
which has traditionally been used to establish the validity of sum-rule techniques. Finally in Section |J wc study 
the phenomenologically challenging case of the quark scalar mesons probed through the non-strange fin current. 
The scalar sector requires an extension of the analysis techniques to a two-resonance phenomenological model, and 
development of criteria to determine when such an extension is necessary. This sensitivity of the Gaussian sum-rule 
to excited states is one of the unique features of the Gaussian sum-rule compared with Laplace sum-rules. The 
phenomenological implications of the results of the Gaussian sum-rule analysis of the fin quark scalar mesons will 
be presented at the end of Section 0. 



2 Conceptual Foundations of Gaussian Sum-Rules 

Consider the basic form of the Gaussian sum-rule (|l|). Leaving aside the method for calculating G(s, r), consider 
the r — > limit of the sum-rule. Using the identity 

= S(t-s) (4) 




we see that 



lim G (s, t) — — p (s) , s > t a (5) 



7T 



Since p(t) is related to hadronic physics quantities, the r — * limit of the Gaussian sum-rule would in principle 
allow direct extraction of hadronic physics from the s dependence of G (s, r) obtained from QCD. 

Although the above limit cannot be achieved in practice, the property (^) demonstrates an important conceptual 
advantage of the Gaussian sum-rule. Consider the equivalent form of the Laplace sum-rule 

oo 

fl(A 2 ) = i [ eX p(-4^]p(t)dt . (6) 



7T J ~' e V A 2 

to 



The A — * limit (again not possible to achieve in practice) would in principle only emphasize the t = region of 
the integration region, and would not reveal the energy dependence of pit). 

Now consider a more realistic case of a Gaussian sum-rule with r > 0. The Gaussian kernel is peaked at t = s, 
and has a width of ^/2t. Thus G (s, r) represents a smearing of p{t) in the region s — ^/2t < t < s + y/2r, centered 
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at t = s. It is therefore reasonable to hope that with r fixed at a reasonable physical value, the s dependence of 
G (s, r) will reveal the rough structure of p{t). For example, if t = m 2 . corresponds to a sharp resonance peak of 
pit), then one would also expect a peak in G (s, r) near s = m 2 .. These expectations are upheld by the analysis of 
the p and scalar quark mesons in subsequent sections. 

The concept of quark-hadron duality, where averaged hadronic quantities are equivalent to the QCD prediction, 
is explicitly manifested in the parameter r which controls the size of the region near t = s over which hadronic 
physics is averaged. As will be discussed below, r is directly related to the energy scale for running quantities in 
the renormalization-group improvement of the theoretical prediction G (s, r), and no constraint on the parameter 
s emerges. The limitations of QCD originating from the renormalization scale then provide a natural duality 
interval for reasonable agreement between hadronic physics and QCD. Finally, we shall see that nonperturbative 
contributions to G (s, r) are suppressed with increasing s, explicitly reinforcing the importance of non-perturbative 
effects in the low-energy region. 

The above aspects of the relation between QCD and hadronic physics illuminated by Gaussian sum-rules are 
obscured in the Laplace sum- rules. The reason for this contrast is that the Laplace kernel exp (— i/A 2 ) exponentially 
suppresses the entire energy region of p(t), while the Gaussian kernel only damps the energy region of p(t) away 
from t = s. Thus the two scales s and r in the Gaussian sum-rule provide a more detailed probe of hadronic physics 
than the Laplace sum-rule. 



2.1 Single-Resonance Analysis of Gaussian Sum- Rules 

Extraction of hadronic properties from the theoretical prediction G (s, r) requires a phenomenological model for 
p(t). The "resonance plus continuum" model is used to represent the hadronic physics phenomenology contained in 
the integral of p(t). In this model, hadronic physics is (locally) dual to the theoretical QCD prediction for energies 
above the continuum threshold t = s 

p(t) = 9 (s - t) p had (t) + 9(t- s ) lmIl QCD (t) (7) 

where sq is the continuum threshold above which hadronic physics and QCD are locally dual. The continuum 
contribution to the integral of (|l|) is denoted by 

°° / 2 \ 

G^( S) r,s ) = I J ^exp f-^-J ImUQ CD (t)dt (8) 

so ^ ' 

Since the continuum contribution is determined by QCD, it is combined with the theoretical quantity G (s,r) 

G QCD (J, r, So ) = G (J, r) - G cont (J, r, s ) (9) 
resulting in a Gaussian sum-rule relating QCD to hadronic physics phenomenology. 

G^ CD (J, r,s ) = iy -j= exp -Z^p- P had (t) dt (10) 

to ^ ' 

Consider the single narrow resonance model for p had (t) widely (but not exclusively) employed in the analysis of 
Laplace sum- rules 

\p had {t) = f r 6(t-ml) (11) 



where m r denotes the resonance mass, and / 2 is the integrated resonance strength. Use of this model in ( |io|) leads 
to a Gaussian sum-rule which relates QCD to the properties of the single resonance. 



GQ° D (s,r, So ) = ^exp 



(s — ml) 2 
17 



(12) 
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The continuum threshold so is constrained by the finite-energy sum rule (|^) , which in the single narrow resonance 
model becomes 



F (s ) = f~ 



(13) 



Since the finite-energy sum-rule must be satisfied for the Gaussian sum-rule to evolve asymptotically to the QCD 
prediction through the diffusion equation (@) (heat-evolution test Q ) , the normalization of the sum-rule (||) is fixed 
by the finite-energy sum-rule constraint, as can be verified by taking the s integral of (fL2|). 



DC 

J G QCD { 



s, t, s ) ds = fl 



(14) 



Thus since the quantity is constrained by the finite-energy sum-rule (^|), the normalization of the quantity 
qQCD ^ T ^ s ^ j g a i reac [y determined. Consequently, the information in the sum-rule ( |l2| ) remaining after imposing 
the heat-evolution test p] is contained in the quantity N^ CD (s, t, sq) normalized to unit area 



N* cd (s,t, So ) 



GQ CD (s,t,s ) 



M (t, s ) 
M (t,s ) = / G QCD (s,T,s Q )ds , 



(15) 
(16) 



leading to a sum-rule independent of the resonance-strength and finite-energy sum-rule. 



N* cd (s,t,s ) 



I 4ttt 



exp 



4r 



2 1 



(17) 



Viewed as a function of s, the phenomenological (right-hand) side of (|17|) has a peak (maximum) which occurs 
at s = m 2 r independent of t, and the peak height 1/v 'Attt at this maximum shows r dependence. These properties 
will be used to find the optimum value of so from which m 2 can be extracted. This optimum value of so can then 
be used to determine the total resonance strength through the finite-energy sum-rule (13). 

The analysis proceeds by fixing r, and finding the value s pea k (t, so) at which N® CD (s, t, so) has a maximum^] 



= 



S = Speafc(l",S ) 

and then determining the peak height 

N P eak (r, s ) = N QCD (s peak (t, s ) , r, s ) 



(18) 



(19) 



The optimum value of So is constrained by the requirements following from the properties of the phenomenological 
side of @: 

1- Speak (t, so) — to 2 , is constant (r-independent), 
2. V A-kt Npeak (t, sq) — 1 (independent of r). 

Thus the optimum value of so, and the corresponding prediction of the resonance mass m 2 , can be determined by 
minimizing a x 2 measure related to the above properties 



(-0) = £ 



Speak -^o) 



l) ( N peak {n, So) \/47rn-l) 



(20) 



1 Our numerical analysis uses the Golden search algorithm JH| to determine s pea fc (r, so) 
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where Tj are equally-spaced points in the r region of interest, and m 2 . is implicitly a function of Sq 



i=i 



^peak 



(n, s ) 



(21) 



The value of so which minimizes the x 2 thus provides the best agreement with properties of the phenomenological 
model and leads to a prediction of the resonance mass m 2 in the single narrow-resonance model. 

Before proceeding with the application of these \ 2 techniques to the p meson, extensions of the narrow resonance 
model will be discussed to determine the magnitude of effects that could arise from resonance widths T and resonance 
shapes. Intuitively, one would expect that if the Gaussian width is much wider than the resonance width, then the 
width effects will be negligible. 

The simplest extension of the unit-area narrow resonance S (t — m~) is the unit area square pulse previously 
used for studying width effects in the Laplace sum-rules H 



1 



P sp (t) 



1 



[6(t- 



<>R 



r T)-6(t-m 



m 



•r)] 



(22) 



which has the following Gaussian image 

ml+m r r 



G sp (s, t) 



where 



1 



1 



4m r r 



exp 



— m r T 
2 



(t-sy 



s — m R + m r T 

V? 



4r 



dt 



(23) 



erf ( § - m l~ m rT 



erf(x) 



dy 



(24) 



Such a resonance model clearly overestimates the effect of resonance widths compared with a Brcit-Wigner which 
is more concentrated about the resonance peak. The expression (|23| ) still has a peak position at s = m 2 as in the 
narrow resonance model, but its peak height is altered to 



G sp 



1 



2m r T 



-erf 



/ m r r\ 


1 




V / 47TT 



1 



m 2 r T 2 
12t 



or 4 



(25) 



As anticipated, the resonance width effects in ( p5| ) diminish with increasing Gaussian width r. Using the full 
expression in terms of the error function, the deviation from the narrow width limit is less than 5% for m r T / < 
0.4, a condition which is easily satisfied for even a wide resonance such as the /o(1370) in the r > 1 GeV 4 range. 
Such uncertainties are well below those associated with the theoretical prediction of G (s, t). For the p meson, the 
resonance width is a completely negligible 0.2% effect even for r = 0.5 GeV 4 . Thus we conclude that for reasonable 
ranges of the the Gaussian width t, the Gaussian sum-rule is insensitive to the effect of resonance widths, and so 
the narrow width phenomenological model is in fact an accurate description of non-zero width resonances in the 
Gaussian sum-rules. 

2.2 Formulation of Gaussian Sum-Rules 

In this section we briefly summarize salient features of the formulation of Gaussian sum-rules presented in 
Consider a dispersion relation with one subtraction constant as needed for the correlation function of vector currents 
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used to probe the p meson]^. 

oo 

n(g 2 ) =n(o)-^/dt 



*(* + Q 2 ) 



(26) 



The difference between the dispersion relation at Q 2 = — s + iA and Q 2 = —s — iA cancels the dependence on the 
subtraction constant n(0). 



n (-s - iA) - II (-§ + iA) _1 . , 
2iA " ~ 1 



p(t) 



(t-iy + A 2 



(27) 



For large A the left-hand side of (j27y is determined by QCD since it evaluates the correlation function well away 
from the physical cut. As with Laplace sum-rules Jl|, the Gaussian sum-rule can be obtained from the dispersion 
relation (J27|) through the Borel transform operator B 



D 



lim 

A 2 /N=4t 



-A 2 



T(N) \dA 2 



N 



(28) 



which has the following property relevant to construction of the Gaussian sum-rule 
1 



B 



A 2 



^ exp (-£) 



Defining the theoretically-determined quantity 

"II(-s-iA) -II(-sH-iA) 



G(s,r) = J r ~B 



iA 



(29) 



(30) 



leads to the Gaussian sum-rule (P after application of B to (p7|). 

An alternative to the direct calculation of the Gaussian sum-rules through the definition of B in (§8|) is obtained 
through an identity relating the Borel and Laplace transform [0 . 



oo 

/ (A 2 ) = J F{T)e~^ T dT = C [F(T)] 
o 

6+ioc 

£" 1 [/(A 2 )]= 5 L / /( A 2 )e^dA 2 



^B[f(A 2 )]=F(T) = C-'[f(A 2 )] , T = (31) 



(32) 



where the real parameter b in the definition (^) of the inverse Laplace transform must be chosen so that / (A 2 ) 
is analytic to the right of the contour of integration in the complex A 2 plane. Using the result (|n|), the Gaussian 
sum- rule (30) can be obtained from an inverse Laplace transform 



G(s,t) 



40FF 



II(-s-iA) -II(-s + iA) 



iA 



1 1 

A\/tFF lixi 



II(-s-iA) -II(-s + iA) 
iA 



exp (A 2 T) dA 2 



(33) 



b-i 



2 Extension to correlation functions requiring further subtraction constants leads to a final result identical to ( pc| ) obtained for a 
single subtraction constant. 
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Respectively mapping the individual terms in (|3^) containing II (—3 ± iA) from the A 2 to the w 
plane results in the following expression for the Gaussian sum-rule 



-s±iA complex 



G(s,r) 



1 1 





(w + s) 2 
4r 


J H(w) exp 


+r 2 





dw 



(34) 



where the contours Ti and T2 are illustrated in Figure |I| 



3 Prediction of the p Mass from Gaussian Sum-Rules 

The p meson is probed through the vector-isovector current correlation function. 

2 = -q 2 > 



n, (Q 2 ) [q»q v - q 2 g^] = i j d 4 x e^ x (O \T [J»(x) J"(0)]| O) , 
J»(x) = - [u(x)^u(x) - d(x)^d(x)] 



(35) 
(36) 



The field-theoretical (QCD) calculation of U v (Q 2 ) consists of perturbative (logarithmic) corrections and QCD 
vacuum effects of infinite correlation length parametrized by the power-law corrections from the QCD vacuum 
condensates M 



u v (q 2 ) = np v ert (q 2 ) + ni ond (q 2 



(37) 



To two-loop order in the chiral limit m u = ma = the perturbative contribution^] in the MS scheme for three active 
flavours is M 



1 



n: ri (Q 2 ) = -^iog 



Q' 1 



1 



(38) 



where v is the renormalization scale. To lowest order, the QCD condensate contributions up to dimension 8 arc 
given by || 



v ^ ' 87T 2 



1 



1 



1 



(C4O4) = — (aG 2 ) - 8n 2 m (qq) 



(39) 
(40) 
(41) 



where SU(2) symmetry and the vacuum saturation hypothesis have been employed. For brevity, we refer to the 
literature j8j for the expressions for the dimension eight condensates, and simply use ( [39] ) to establish a normalization 
consistent with ||. 

To evaluate the Gaussian sum-rule consider the contour C(R) in Figure |2|. The quantity Tl(w) is analytic within 
and on C(R); consequently 



1 1 

/47TT 27TI 



<j> H(w) exp 



C(R) 



(w + §y 

4r 



dw 



(42) 



In the limit as R — > 00, the integrals along T3, and T$ all approach zero for the QCD expressions given in 
(38 1 |39| ). Furthermore ^i(R) and fi(-R) respectively approach Ti and T2 as R — * 00. We therefore obtain the 

3 A divergent constant has been ignored since like the subtraction constant 11(0), it vanishes in the formation of the Gaussian sum 



rule. 
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following expression for the Gaussian sum-rule 
1 1 



+ 



lim 

/47TT 27TI R— >oo 



H(w) exp 



r c +r e 



(w + §y 

4r 



dw 



1 



1 



/4ttt 2iri 



[n (ie i7r ) - II (ie'™)] exp 



_L J_ / ee ^n ( ee ie ) exp - 

— 7T * 



(ee* e + s) 
47 



4r 



(43) 



The perturbative contributions from the T e contour (8 integral) in (|43) will be zero in the e — > limit, leav- 
ing only the integral of the discontinuity across the branch cut [i.e. JmW (t)] to determine the perturbative 
contributions to the Gaussian sum-rule 



G pert ^ T j 



1 



'47TT 



1 

8^ 



a(y) 



exp 



4r 



I 



16tt 2 



a(i/) 



1 + erf 



2^ 



, (44) 



Since the current d36| ) is renormalization-group invariant, Imil per *(t) satisfies a homogeneous renormalization-group 
equation. The result of renormalization-group improvement for the Gaussian sum-rules can be inferred from the 
general structure of perturbative corrections which take the form 



log" 



t 



exp 



At 



dt = v 2 / dx log™ (x) exp 



2v^ 2^? 



2 

7^ 



(45) 



The functional dependence of the perturbative corrections to the Gaussian sum-rules expressed by the function 
H demonstrates that v scales with y/r. Thus, the solution of the renormalization-group equation satisfied by the 
perturbative contributions is the replacement of a(v) with the running coupling constant and identifies the scale 
v 2 = t/t. This implies the existence of low-energy boundary (lower bound) on r, but places no restriction on the 
scale s. 

The QCD condensate contributions ( |39| ) to II (Q 2 ) do not have a branch discontinuity, so their contribution to 
the Gaussian sum-rule arises solely from the contour T e (6 integral) in p3), and can be evaluated using the result 



2iri J w 1 



exp 



(w + s) 
4t 



dw = lim 



™-»o (n - 1)! dw™- 1 



exp 



(w + §y 

4r 



n = 1,2,3. 



(46) 



Expressions (|44j) , (46), d39), and (|43j) lead to the Gaussian sum-rule of vector currents. 



G(s,r) 



1 



16tt 2 
+ 



a(y/r) 



1 



64n 2 T\/Trr 



1 + erf 

'2 



2v^ 



327r 2 ry / 7rr 



exp 



s 

"4t 



2r 



1+ exp (CJC7 6/ 



4r 



1287r 2 r 2 ,/^F 



-l + -]exp 



■5; ) (CM 



(47) 



Agreement between (|47|) and |5| provides a useful check on the conventions established in (|34|). We also note that 
the QCD condensate contributions to p^ ) can be recovered by expanding the Gaussian kernel in (|I]) about t = 
(or alternatively in a series about any other value of t, or in a series of Hermite polynomials), and then using 
the relation between the QCD condensates and the FESR family ^ (in addition, see equation (6.15) in reference 
||). This demonstrates that in principle (e.g. knowledge of higher-dimension condensate contributions), sum-rules 
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based on the Gaussian (or Laplace) kernel contain the same information as the whole FESR family. However, 
as we discussed previously and will be demonstrated below, the Gaussian kernel arranges this information more 
effectively, enhancing the predictive power compared with other sum-rule kernels. 

The QCD continuum contributions (ra) arising from the perturbative corrections to the vector correlator are 



G cont (s,r, So ) 



1 



16?r 2 



1 



a(Vr)\ 



1-erf 



so 



2^ 



and hence the theoretically determined Gaussian sum-rule quantity (m for the vector currents is 



G QCD { 



1 



S,T, S J 



16tt 2 
+ 



1 



erf 
s 2 \ 



( 



6Att 2 t^/ttt 
s 



2r J 



exp 



erf 



so 



s 

"It 



s 

6r 



exp 



s 

"It 



2V^ 

(ClOl) 



327T 2 1 



exp 



s 

"47 



(48) 



(49) 



Finally, the quantity Mo (r, so) in (|l6|) required for the normalized vector-current Gaussian sum-rule (JTJ^) is given 

by 



M (t, Sq )= J G® CD (s,T, So )ds=-^(l 



a(y/r) 



(50) 



The non-perturbative QCD condensate contributions in (|49|) are exponentially suppressed for large s. Since s 
represents the location of the Gaussian peak on the phenomenological side of the sum-rule, the non-perturbative 
corrections are most important in the low-energy region, as anticipated by the role of QCD condensates in relation 
to the vacuum properties of QCD. This explicit low-energy role of the QCD condensates clearly exhibited for the 
Gaussian sum-rules is obscured in the Laplace sum-rule because the peak of the Laplace exponential kernel in ([|) 
is always located at t = 0. 

Several QCD parameters must be specified before carrying out the phenomenological analysis of the Gaussian 
sum-rule of vector currents. The running coupling for three active flavours to two-loop order is 



a, {v 2 ) = J_ 0i \ogL 
7T l3 a L (p oL f 



L = log 



V 

A 2 



0i 



(51) 



with A -jjg sa 300 MeV for three active flavours, consistent with current estimates of a s (M T ) Jl0|, |TT| and matching 
conditions through the charm threshold |L3| . For the gluon condensate we use the central value determined in M 



(aG 2 ) = (0.045 ± 0.014) GeV 4 , 
and PCAC jia] is used for the quark condensate. 



, (qq) = ^m u (uu) + ^m d (dd) = -~f%ml , f n = 93 MeV 



(52) 



(53) 



Deviations of the dimension-six condensate from the vacuum saturation value are parameterized by the quantity 
f vs which could be as large as f vs — 2 P, [l4[ 



, 896 3 ., .- a 



896 



(COT = -f vs ^^a {{qq)Y = -f vs ^ir a (1.8 x 10- 4 GeV b ) , f vs 

Finally, we use the reference |9j value for the dimension-eight condensate. 
(ClOl) = (0.40 ±0.16) GeV 8 



1.5±0.5 



(54) 



(55) 
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The x minimization techniques outlined in Section 2.1 [see equations (|20|) and (21)] are now employed to 
determine the optimum continuum threshold sq and the corresponding prediction of the p meson mass m p , demon- 
strating the ability of Gaussian sum-rules to predict resonance properties. Using equally-spaced points in the 
range^J 0.5 GeV 4 < r < 4.0 GeV 4 for determining % 2 (so) in (|20|), and with inclusion of the uncertainties of the 
QCD condensates given in (p2]-p35|) we obtain the following results for m p and Sq 



rn, 



(0.75 ± 0.07) GeV , s = (1.2 ± 0.2) GeV 2 (56) 



in excellent agreement with the measured mass m p = 770 MeV. 

A detailed examination of the required equality (0) between the Gaussian sum-rule and the phenomenological 
single resonance model can be obtained by using the x 2 -predicted values of m p and so as input into the single- 
resonance plus continuum model, and examining the s, r dependence of the sum-rule in comparison with the 
phenomenological model. Figure |^ compares this s dependence of the sum-rule and phenomenological model for 
selected t values in the region used to define the \ 2 . As would be anticipated by its construction, the \ 2 optimization 
procedure leads to excellent correspondence between theory and phenomenology for the s pea k position and peak 
height. However, the astounding agreement between the s dependence of the sum-rule and phenomenological model 
provides strong evidence for the ability of Gaussian sum- rules to predict hadronic properties. In the next section, 
Gaussian sum-rules will be employed to study the (non-strange) nn scalar mesons, a more challenging channel for 
hadronic physics phenomenology. 



4 Gaussian Sum-Rule Analysis of the Quark Scalar Mesons 

The interpretation of the scalar mesons is a challenging problem in hadronic physics since a variety of interpretations 
exist for the lowest-lying scalar resonances (a or / (400_- 1200), / (980), /o(1370), /o(1500), a (980), a (1450) @) 
including conventional quark-antiquark (qq) states, KK molecules, gluonium, four-quark models, and dynamically 
generated thresholds [l5| [l6| 

The relevant correlation function for the I = 0, 1 non-strange fin scalar mesons is 

n s (Q 2 ) =i [ d A xe^ x (0\T[J s (x)J s (0)}\0) , Q 2 = -q 2 > (57) 



where, in the SU(2) limit \m q = (m u + md)/2] for isoscalar (/ = 0) and isovector (/ = 1) currents, 

J s (x) = m q [u(x)u(x) + d(x)d{x)] /2. (58) 

We note that the factor of the quark mass is necessary for a renormalization-group invariant current. Correlation 
functions of scalar and pseudoscalar currents are unique since they receive significant contributions from instantons 
|T^ , |l9| . In contrast to QCD condensates which represent vacuum effects of infinite correlation length, instanton 
contributions represent finite correlation-length QCD vacuum effects, and are the only known theoretical mechanism 
that distinguishes between the 1 = and 1=1 correlation functions in the presence of SU(2) flavour symmetry 
]2C| . Thus the theoretical calculation of the scalar current correlation function consists of perturbative, condensate, 
and instanton contributions 

n s (q 2 ) = nf rt (q 2 ) + n c s ond (q 2 ) + ir s nst (Q 2 ) . (59) 

At two-loop order and to leading order in the quark mass for three active flavours in the MS scheme, the perturbative 
contributions in ( |59| ) are pl|]p| 



. . „. 3m 2 Q 2 /Q 2 



a f , Q z \ 17 



7 " log + T 



(60) 



4 Altering the range to 1.0 GeV 4 < r < 4.0 GeV 4 has minimal effect on the predictions. 

5 A field-theoretical divergence proportional to Q 2 is ignored since it vanishes after application of B in the formation of the Gaussian 
sum-rule. 



10 



To leading order in m q and a, the QCD condensate contributions in ( p9| ) are |2^, [23] [[]: 



ncond //~)2 \ 



Q 2 



(m q qq) 



16tt 



<a,G 2 ) 



(<7|0f> = 7T« S 



(u7„A u + d 7M A a d) J] q^X a q 



u.d.s 



(61) 
(62) 

(63) 



As for the vector current correlation function, the vacuum saturation hypothesis jlj] in the SU(2) limit 
(dd) = (qq) provides a reference value for 



(c s 6 o 6/ 



fvs 2j a ' 



(m 2 ) 



-f vs 1.8 x lCT 4 GeV b 



(64) 



As mentioned earlier, the perturbative and condensate contributions to the scalar correlation function are indepen- 
dent of J, and hence do not distinguish between the isoscalar and isovector channels. 

Using symmetry properties relating the instanton contributions to pseudoscalar and scalar correlation functions 
PQl combined with the instanton liquid model results for the pseudoscalar correlation function |l], |l9| , the instanton 
contributions to the scalar correlation function are 



nr* (Q 2 ) = (-i) J 3m '° 2 



Air 2 



(65) 



where K\ is a modified Bessel function J2j| , and p c — 1/ (600 MeV) is the (uniform) instanton size in the instanton 
liquid model [[l9| . The explicit factor depending on / in ( |65| ) is the only theoretical contribution that distinguishes 
between the isovector and isoscalar channels. 

Calculation of the Gaussian sum-rule for the scalar currents proceeds as in Section Ej. Using the following 
asymptotic property of the modified Bessel function (2j| 



K x (z) 



\z\ > 1 , |arg(z)| < 



3tt 



(66) 



one can verify that the QCD expressions given in (|6(J, £35|) uphold the transition from ([12]) to ( fl3| ) developed 
for the vector currents. Thus we can calculate the Gaussian sum-rule for scalar currents from (B3[). 



As with the vector channel, the perturbative contributions (60) from the contour T e [0 integral) in (M3J) will 
be zero in the e — > limit, leaving only the integral of the discontinuity across the branch cut [i.e. lmn? e {t)\ to 
determine the perturbative contributions to the Gaussian sum-rule. 



1 3to 2 (Jt) 
G pert (s,T) = 1 X 



/AlTT 167T 2 



t 1 



17 a (7?) 



- 2 



*l0g -F= 



exp 



At 



dt 



(67) 



This expression has already been renormalization-group improved by the replacement of v 2 = yfr and implicit 
identification of a and m q as running quantities at this same scale. The running coupling has already been given 
in (|5l|), and the running quark mass for three active flavours at two- loop order is 



290 1 256 log L 
729 L ~ 729 L 



lOB 



A 2 J 



(68) 



where rh q is the renormalization-group invariant quark mass parameter. 

As for the vector channel, the QCD condensate contributions (|6l]) do not have a branch discontinuity, so their 
contribution to the Gaussian sum-rule arises solely from the contour T e (6 integral) in (Q), and can be evaluated 
using (46). 
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The instanton contributions to the Gaussian sum-rule follow from (J33|) 



■ 4 i 3m 2 . 1 1 

G mst (s,r) =(-l) 7 q 



4-7T 2 J Attt 2iri 



j 3m 2 1 1 



+ (-!)' -« 



47T 2 J4t7t2it 



eV 29 



2\ / (t-sy 

exp 



At 



(p c \/ee 



J8 



exp 



4r 



d6> 



(69) 



Simplification of (^39|) requires the following properties of the modified Bessel function K\(z) [^|| 
Ki(z)~- , z^O 



Ki (z) = 



-f H[ 1] (ze"/ 2 ) , -7T < arg(z) < f 
-f-ff x (2) [ze-^l 2 ) , -f < arg(z) < tt 



(70) 
(71) 



where H[^(z) — J\{z) + zYi(z) and h[ 2 \z) = Ji(z) — iYi(z). The asymptotic behaviour ([70[) implies that the 6* 



integral of (66) will be zero in the e — > limit and the identity ( |7l|) allows evaluation of the discontinuity in the t 
integral of (69), leading to the following instanton contribution to the Gaussian sum-rule 



r 3m 2 1 

G mst (s,r) = -(-!)' g 



57T JAtTT 



exp 



At 



dt 



(72) 



Expressions (|67j), (46), (|6l]), ( |43| ) and ( |72| ) lead to the Gaussian sum-rule of scalar currents 



\jAlTT 16 
-(-1) 



exp 



(t-s) 2 



At 



dt 



! 3m 2 1 



8ir \J Attt 

2 \ 



/ t Ji (p c Vt) Yi (p c \/<) exp I - 
o v 

2\/ttt Atwttt 



At 



dt 



(73) 



Comparison of ( |72j ) and (||) indicates the existence of an instanton continuum contribution 

TT OTT 



(74) 



Combined with the perturbative continuum devolving from ( |60| ) we obtain the continuum contribution to the 
Gaussian sum-rule of scalar currents 



G cont (s, t) 



1 3m 2 (V?) 
/Attt 167T 2 



3m 2 1 



til 



3 TT 



Hog 



exp 



(t-sf 

At 



dt 



57T J Attt 



tJi (pcVi) Yi (pcVi) exp - 



At 



(75) 



dt 
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and hence the theoretically determined Gaussian sum-rule (0) for the scalar currents 



c«c» (4 , r) 7 r t (, + _ 2 f^« ,„ g ( « 

y/Airr 16tH jL\ 3 tt / tt \^/t 

n 



exp 



(*-3) 2 ' 
4r 



dt 



M) 



/ 3m 2 i 



: J tJi (pcVt) Yi (pcVt) exp ( -^: s) 

n \ 



dt 



(76) 



4rx/7rr 



Finally, the quantity Mo (r, so) in 

by 



( |l6| ) required for the normalized scalar-current Gaussian sum-rule ( |l5| ) is given 

Mo (t ^ (l + f ^ | _ | (2 , og ( * ) _ l} 

3 2 S ° 
n 



(77) 



All theoretical contributions in ([76]) and ( |77j ) are proportional to m 2 , and hence proportional to the m 2 through 
(68). Consequently, the normalized Gaussian sum-rule N® CD will be independent of m 2 which is clearly advanta- 



geous given the uncertainty in determinations of the light-quark masses [jlOj . 

If the Gaussian sum-rules for scalar currents are analyzed using the single-resonance ^-optimization techniques 
of Section 2.1 we find poor agreement between the theoretical prediction and the single resonance phenomenological 
model for the optimized values of the mass and continuum. Figures || and || show that the peak height of the single 
resonance model is consistently larger than the theoretical prediction, and that the peak position of the theoretical 
prediction shows some r dependence. Furthermore a discrepancy is seen at the tails of the distributions where 
the single-resonance model becomes smaller than the theoretical prediction. This latter point indicates that the 
width of the theoretical distribution is larger than that of the single resonance model. 

A quantitative measure of this behaviour can be found in higher-order moments of the Gaussian sum-rule. 



M k (T,s )= J s k G QGD (s,T,s )ds , k = 0,1,2,3... 

— OO 
OO 

Mfc ( r ' s °) = J s k N Q CD ( Sl T lSo )ds , k = 0,1,2,3... . 



M (r, s ) 



(78) 
(79) 



Of particular importance are the moments related to the quantities a 2 and A2 defining the width and asymmetry 
of the distributions 



a 2 = 



M2 



Mi 
Mo 



_ M 3 M 2 Mi / mi 
A2= Mo Mo Mo + ! ~~ 



Mi 
Mo 



(80) 
(81) 



6 We note that peak drift is absent in the equation (^) expression for width effects, and furthermore unrealistically largi 
widths would be required to accommodate the observed discrepancy at the peak height. 



c resonance 
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For the scalar channel, the relevant higher-order moments for our analysis are given by 



M 1 (r, s ) 



3™ 2 (Vr) 



16tt 2 



1 + 17 a(V?) \ «| _ 3r^(^)q(y7) 2s 3 , / /_sq. 

3 7T / 3 167T 2 7T 9 \ VV^T" 



- 1 



!)' ^ / * 2 ^1 (poVi) Y X (f> c St) dt - TO 2 {C s s 



M 2 {t, s ) 



3™ 2 (^) / 17a(Vr) 



16tt 2 

.,2 



3™ 2 a (V?) 



16tt 2 
3m 2 



7T 

so 



•so 



s^ 4 log -1 +8r 21og -1 



so 



M 3 (r, s ) 



o 

3m 2 (V?) 



16tt 2 

.,2 



1 17 a(V?) 

3 7T 



3m|(Vr)a(V^ 3 

16^ 2 7T ° 

Sn 

3m 



|(51og(il)-l) + ±r(:no,(-^)-, 



so 





For completeness, we also give the moments needed to calculate tr 2 for the vector channel. 



M a (r, s ) 



8tt 2 



7T / 



8tt 



Consider a two-resonance extension of the narrow resonance model (ill]) so that (fL2) becomes 



fi 



'Attt 



exp 



4r 



f'2 



'Attt 



exp 



4r 



The corresponding expression for the normalized Gaussian sum-rule then becomes 



\>4lTT 



r\ exp 



4t 



+ r 2 exp 



4r 



where 



./I 



r l = /2 + /| . r2 = r/Z r.i • ''I - r - = 1 



/l 2 + /I 

If we calculate moments of (|87|), we find that in the two-resonance model 

M l 2 2 1 / 

— = r 1 m l + r 2 m 2 = - (z + ry) 

2 o / 2 2\ 2 1 2 /i 2\ 

cr - 2r = rir 2 (rrij - m 2 J = -y (1 - r J 

3 1 -j 

M = r x r 2 (r 2 - n) (m 2 - m 2 ) = -^H/ 3 (l - r 2 ) 
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where 



n - r 2 



2 i 2 

z = m 1 + m 2 



(93) 
(94) 
(95) 



We see that in the single resonance limit, a 2 = 2t and A 2 = 0, and hence a clear signature of the existence of two 
resonances is er 2 — 2r ^ 0. f\ 

In Figures || and 0, a 2 is plotted as a function of r for the optimum values of sq resulting from the single- 
resonance x 2 analyses of the 7 = 0,1 scalar currents. From Figures ^ and fj] it is evident that a 2 — 2r is significantly 
different from zero, indicating the presence of a second resonance in both the I = 0, 1 scalar channels. In contrast, 
Figure H for the vector channel shows that a 2 differs insignificantly from 2t. Consequently, there is no evidence that 
a second vector resonance has enough strength to necessitate extension of the phenomenological model, implying 
that any such states are weak enough to be absorbed into the continuum, and providing a validation of the single- 
resonance analysis. Compared to the Laplace sum-rules which exponentially suppress excited states, the potential 
sensitivity to excited states is a unique feature of the Gaussian sum-rules. 

How might we extract predictions of the hadronic parameters {mi, r%, m 2 , r 2 } (or equivalently {r, y, z}) from the 
Gaussian sum-rules (|l^) in a two narrow resonance model? The key to solving this problem lies in the behaviour of 
the peak position s pea k(''") So) of (s, t, sq). As noted in Section 2.1, in the single narrow resonance model, the 

phenomenological side of (|l7]) possesses a r-independent peak located at s = m 2 . In the double narrow resonance 
model (p8|), the situation is more complicated. For any fixed value of r sufficiently large compared to m 2 — m\ (in 
practice, this is indeed the case), the right-hand-side of (|8^) exhibits a single peak; however, the position of this 
peak varies with r, and it is this r-dependent peak drift which enables us to extract hadronic parameters in the 
double narrow resonance model. 

To determine an expression for the phenomenological peak drift, we find the location of the peak by differen- 
tiating the right-hand-side of ( |88| ) and setting the result equal to zero. In terms of the parameters (p3|-p5|), this 
yields 



(r+l)(i 



(r-l)(* 



exp 



y{z-2s) 



4r 



= 0. 



(96) 



Unfortunately, we cannot explicitly solve (J9 1 
series expansion in the parameter 1/r: 



A B C 

A+- + -, 

T 



for s as a function of r, so we instead compute the peak drift as a 



(97) 



where A, B, and C are functions of the parameters {r, y, z}, and the quantities {B, C} will be non-zero except in 
the single resonance limit (y = 0) or in the case of equal resonance strengths (r = 0). Explicit numerical testing in 
worst case scenarios involving two resonances with masses in the 1-2 GeV range show that the first omitted term 
in ( |97j ) [i.e. D/t 3 } is negligible provided that we restrict our attention to r > 2 GeV 4 . 

Wc define the following \ 2 function which measures quantitatively the discrepency between the theoretical peak 
drift contained in 7V^ CD (s, t, sq) and the phenomenological peak drift approximated by ( p7| ) 



N 



2 (*o) = £ 



3peafc(so, U) - A(s ) 



B{s ) C{s ) 



-i 2 



n 



(98) 



where A(sq), B(sq), and C(sq) are obtained by solving the linear system of equations defined by the ^-minimizing 
conditions 



dx 2 = dx^ _ 
dA~ dB ~ dC 



0. 



(99) 



7 Note that A2 = for resonances of equal strength r\ = r2, and hence A2 = does not necessarily imply a single resonance scenario. 
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Minimizing ( |98| ) numerically using the Golden search algorithm || determines an optimum value for sq- At this 
optimum so, the moments ([78]) exhibit negligible dependence on r as explicitly occurs in the double resonance 
model, and thus the moments are well approximated on the entire r region of interest by averaged values 

Mfco) = I ' Mi(r, s Q )dT. (100) 



Next, we calculate a 2 — 2t and A 2 by substituting the average moments (100) into (79) and (80). Finally, we 



substitute these quantities along with Mo and Mi into the following inversion formulae which follow easily from 



(90-32 



2— + ^-^ — (101) 
M cj 2 -2t V ' 



r«3 (102) 

(T z — It 
M 

r = , (103) 
^JA\ + 4(ct 2 - 2r)3 

to arrive at our predictions for the hadronic parameters {r, y, z} in the double narrow resonance model. 

These techniques for analyzing the Gaussian sum-rule for a two-resonance phenomenological model can now be 
applied to the fin quark scalar mesons. The range 2.0 GeV 4 < t < 4.0 GeV 4 will be used to define the \ 2 ( p8|) since 
it leads to acceptable convergence of the series solution ( |97|) for the peak drift as well as minimizing uncertainties 
from the (large) higher order perturbative coefficients in (|7(j). Table [l] shows the Gaussian sum-rule prediction of 
the quark scalar resonance parameters and continuum resulting from the two-resonance ^-optimization procedure. 
Figures [| and |Io| compare the s dependence of the Gaussian sum-rule and the two-resonance model with the Table 
[l] predictions of the resonance parameters. The nearly-perfect agreement between the theoretical prediction and 
two-resonance phenomenological model exhibited in Figures p| and [To| is clearly a vast improvement compared with 
the one-resonance model results displayed in Figures || and||, and provides compelling evidence for the ability 
of Gaussian sum-rules to predict resonance properties. The absence of any indication of disagreement between 
the theoretical prediction and two-resonance model suggests that any further resonances are weak enough to be 
contained in the QCD continuum. 

The prediction of the resonance parameters for the fin scalar resonances are of phenomenological interest. The 
1 = states in Table [j] suggest that the Jo (980) has a significant fin component, with a weaker excited state 
with an fin component lying in the vicinity of the /o(1370) and /o(1500). The absence of a light (substantially 
less than 1 GeV) 1 = state indicates a decoupling of a light a from the fin scalar currents, suggesting a non-fin 
interpretation for the a. The 1=1 states in Table |l| similarly suggest a non-fin interpretation of the ao(980), and 
are consistent with identification of the ao(1450) as the lightest 1=1 state with a significant fin content. Our 
prediction of an / = 1 state near 1.8 GeV suggests identification of the the 1=1 state V(1775) |I(J as a scalar 
state. Finally, we note that the Gaussian sum-rule predictions for the lightest / = 0, 1 fin scalar mesons confirm the 
results obtained from Laplace sum-rules [^7| which, because of the their exponential suppression of excited states, 
are best suited to probing the ground states in these channels. 



/ 


mi 


m 2 


n 


T2 


so 





0.97 GeV 


1.43 GeV 


0.634 


0.366 


2.35 GeV 2 


1 


1.44 GeV 


1.81 GeV 


0.570 


0.430 


3.50 GeV 2 



Table 1: Results of the two-resonance analysis of the Gaussian sum-rules for the I = 0,1 nn quark scalar mesons 
resulting from the x 2 -optimization procedure. Central values of the QCD parameters have been employed. 

As a check on the methods developed for analyzing the two-resonance case, we have also performed a least- 
squares fit between the complete two-dimensional (s, r) dependence of the two-resonance phenomenological model 
and the theoretical prediction. The resonance parameters and continuum parameters obtained from this analysis 
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differ from those of Table]!] by at most 10%, demonstrating that the peak-drift method reproduces a full fitting proce- 
dure with minimal computational expense, and provides unambiguous results compared to a full multi-dimensional 
least-squares fit. At the very least, the peak-drift method can be used to initialize a full multi-dimensional least- 
squares fitting procedure. 



5 Conclusions 

In this paper we have developed techniques for using QCD Gaussian sum-rules to predict hadronic resonance 
properties, and discussed the conceptual advantages of the Gaussian sum-rule in providing a natural duality interval 
for the relation between QCD and hadronic physics. Methods for analyzing a single resonance model based on the s 
peak position and peak height of the theoretical prediction for the Gaussian sum-rule were developed and applied to 
the p meson as a test case. Motivated by the accuracy of the p mass prediction and exceptional agreement between 
the theoretical prediction and single resonance phenomenological model (see Figure ||), the more phenomenologically 
challenging case of the non-strange nn quark scalar mesons was considered. However, the single resonance analysis 
of the Gaussian sum-rule of the quark scalar currents revealed a discrepancy between the theoretical prediction 
and phenomenological model (see Figures | and |). 

This discrepancy between the theoretical prediction and single resonance model for the scalar currents is ad- 



dressed by including a second resonance in the phenomenological model. The quantity cr 2 , defined in (80) by 
moments of the Gaussian sum-rule, provides a criterion for evaluating the necessity of including a second resonance 
in the phenomenological model. Figures |6|-|| illustrate that the single resonance model is validated for the vector 
case, but is insufficient for the scalar currents. Compared with the Laplace sum-rule, the potential sensitivity of 
the Gaussian sum-rule to excited states is a novel feature of the Gaussian sum-rules. 

The methods developed for analyzing the Gaussian sum-rules in a single resonance phenomenological model 
can be extended to the two-resonance model by studying the r dependence of the s peak position (peak drift) of 
the theoretical prediction for the Gaussian sum-rule. When applied to the Gaussian sum-rule of nn quark scalar 
currents, we find the discrepancy between the theoretical prediction and phenomenological model observed in the 
single-resonance analysis is resolved by inclusion of a second resonance, and the resulting agreement between the 
theoretical prediction and phenomenological model (see Figures | and 0) is astounding. 

The phenomenological results for the fin quark scalar mesons (see Table [j]) provide valuable information for 
interpretation of the scalar mesons. In particular, the Gaussian sum-rule analysis is in excellent agreement with 
the identification of the ao(1450) as the lightest quark scalar meson with a significant nn content, supporting 
the conclusions of , and suggests that a light (<S 1 GcV) 1 = scalar meson does not have a significant 

fin component. The 1=1 excited state in Table [l] suggests identification of the A(1775) |l(J as a scalar state 
containing nn. 

In summary, we conclude that the Gaussian sum-rules are not only useful in their well-established connection 
with the finite-energy sum-rule constraint on the continuum threshold [0 , but are a valuable technique for studying 
hadronic physics, and should be considered as a viable and complementary alternative to Laplace sum-rules. 
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Figure 2: Closed contour C{R) used to calculate the Gaussian sum-rule defined by (|34|). The inner circular segment 
r e has a radius of e, and the circular segments T3, T4 and T$ have a radius R. The wavy line on the negative real 
axis denotes the branch cut of II (id), and the linear segments of the contour above and below the branch cut are 
denoted by r c . The contour ri(-R) is that portion of Ti (see Figure [l]) which lies in the interior of a circle of radius 
R centred at — s, and similarly for T2(R). 
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Figure 3: Comparison of the vector-current theoretical prediction for the N® CD (s, t, so) with the single resonance 
phenomenological model using the x 2 -optimized values of the resonance mass m p = 0.773 GeV and continuum 
so = 1.22 GeV 2 . The r values used for the four pairs of curves, from top to bottom in the figure, are respectively 
r = 0.5 GeV 4 , r = 1.0 GeV 4 , r = 2.5 GeV 4 , and r = 4.0 GeV 4 . Note the almost perfect overlap between the 
theoretical prediction and phenomenological model. Central values of the condensate parameters have been used. 
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Figure 4: Comparison of the 1 = scalar current theoretical prediction for N^ CD (§, t, s ) with the single resonance 
phenomenological model using the x 2 -optimized values of the resonance mass m = 1.21 GeV and continuum s = 
2.60 GeV 2 . The r values used for the three pairs of curves, from top to bottom in the figure, are respectively 
t = 2.0 GeV 4 , r = 3.0 GeV 4 , and r = 4.0 GeV 4 . The phenomenological model is consistently larger than the 
theoretical prediction near the peak, but is consistently smaller than the theoretical prediction in the tails. Central 
values of the condensate parameters have been used. 
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Figure 5: Comparison of the 1=1 scalar current theoretical prediction for N^ CD (§, t, s ) with the single resonance 
phenomenological model using the x 2 -optimized values of the resonance mass m — 1.68 GeV and continuum s = 
3.9 GcV 2 . The r values used for the three pairs of curves, from top to bottom in the figure, are respectively 
r = 2.0 GeV 4 , r = 3.0 GeV 4 , and r = 4.0 GcV 4 . The phenomenological model is consistently larger than the 
theoretical prediction near the peak, but is consistently smaller than the theoretical prediction in the tails. Central 
values of the condensate parameters have been used. 
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Figure 6: Plot of a 2 for the theoretical prediction (solid curve) compared with o 2 — 2r for the single-resonance 
model (dashed curve) for the I — scalar channel using the x 2 -optimized value of the continuum s = 2.6 GeV 2 . 
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Figure 7: Plot of a 2 for the theoretical prediction (solid curve) compared with a 2 — 2r for the single-resonance 
model (dashed curve) for the I ~ I scalar channel using the x 2 -optimized value of the continuum s = 3.9 GeV 2 . 
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Figure 9: Comparison of the 1 = scalar-current theoretical prediction N® CD (s,t, s ) with the two-resonance 
phcnomcnological model using the x 2 -optimized values of the resonance masses (mi = 0.97 GeV, m 2 = 1.43 GeV), 
relative resonance strengths (r\ — 0.634, = 0.366) and continuum sq = 2.35 GeV 2 . The r values used for the four 
pairs of curves, from top to bottom in the figure, are respectively r = 2.0 GeV 4 , r = 3.0 GeV 4 , and r = 4.0 GeV 4 . 
Note the almost perfect overlap between the theoretical prediction and phenomenological model. Central values of 
the condensate parameters have been used. 
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Figure 10: Comparison of the I — 1 scalar current theoretical prediction N® CD (s,t, s ) with the two-resonance 
phcnomcnological model using the x 2 -optimized values of the resonance masses (mi = 1.44 GeV, m 2 = 1.81 GcV), 
relative resonance strengths (r\ — 0.570, = 0.430) and continuum sq = 3.50 GeV 2 . The r values used for the four 
pairs of curves, from top to bottom in the figure, are respectively r = 2.0 GeV 4 , r = 3.0 GeV 4 , and r = 4.0 GeV 4 . 
Note the almost perfect overlap between the theoretical prediction and phenomenological model. Central values of 
the condensate parameters have been used. 
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